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Abstract 

We consider the distributions of the lengths of the longest weakly increasing and strongly 
decreasing subsequences in words of length N from an alphabet of k letters. (In the limit as 
k — > oo these become the corresponding distributions for permutations on N letters.) We 
find Toeplitz determinant representations for the exponential generating functions (on N) 
of these distribution functions and show that they are expressible in terms of solutions of 
Painleve V equations. We show further that in the weakly increasing case the generating 
function gives the distribution of the smallest eigenvalue in the k x k Laguerre random matrix 
ensemble and that the distribution itself has, after centering and normalizing, an iV — > oo 
limit which is equal to the distribution function for the largest eigenvalue in the Gaussian 
Unitary Ensemble of k x k hermitian matrices of trace zero. 

I. Introduction 

The last decade has seen a flurry of activity centering around connections between com- 
binatorial probability on the one hand and random matrices and integrable systems on the 
other. From the point of view of probability theory, the quite surprising feature of these 
developments is that the methods came from Toeplitz determinants, integrable differential 
equations of the Painleve type and the closely related Riemann-Hilbert techniques as they 
were applied and refined in random matrix theory. Using these methods new, and appar- 
ently quite universal, limiting laws have been discovered. One of the aims of this paper is 
to make these methods accessible to a wider audience. Our story begins with a theorem of 
Gessel |16| . (There are earlier signs of these connections; see Regev fl30| .) 



Let Sn be the symmetric group on N letters and give each permutation a G Sn proba- 
bility 1/N\. Denote by ^(cr) the length of the longest increasing subsequence in a and 

F P (n; N) = Prob (£ N (a) < n) . 
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Then it is a corollary of Gessel's theorem and the Robinson-Schensted-Knuth correspon- 
denceR that 

u oo 

Y, F p(n;N)—=D n (t), (1.1) 

N=0 iV ' 

where D n (t) is the determinant of the the n x n Toeplitz matrix with the symbol e^^ z+z \ 
(Recall that the i, j entry of a Toeplitz matrix equals the % — j Fourier coefficient of its 
symbol.) 

It is in this work of Gessel, expressing the (exponential) generating function of Fp as a 
Toeplitz determinant, and the subsequent work of Odlyzko et al. |26|] and Rains [pSQ , that 
the methods of random matrix theory first appear in RSK type problems.^ 

Starting with this representation, Baik, Deift and Johansson [0, using the steepest de- 



scent method for Riemann-Hilbert problems [ 121 , derived a delicate asymptotic formula for 



D n (t) which we now describe. Introduce another parameter s and suppose that n and t are 
related by n = [2t 1/2 + s£ 1/6 ]. Then as t — ► oo with s fixed one has 

lim e- t D 2tl/ 2 +st i /(i {t) = F 2 {s). 

t—*oo 

Here F 2 is the distribution function defined by 

F 2 (s) = exp (- [°°(x - s)q(x) 2 dx] (1.2) 



where q is the solution of the Painleve II equation 

q" = sq + 2q 3 



satisfying q(s) ~ Ai(s) as s — > oo^\ Using a dePoissonization lemma due to Johansson [^0 
these asymptotics for D n (t) led Baik, Deift and Johansson to the limiting law 

UmProb /M^_2v^ < \ =FjW . 

It is a remarkable fact that this same distribution function F 2 was first encountered by 
the present authors JR] in random matrix theory where it arises as the limiting law for 



the normalized largest eigenvalue in the Gaussian Unitary Ensemble (GUE) of Hermitian 
matrices. More precisely, we have for this ensemble [p4| 



lim Prob ((A max - V2N)V2N 1/e < s) = F 2 (s). (1.3) 



Here we see a connection with integrable systems — the appearance of a Painleve II func- 
tion. Yet another connection is that D n (t) itself has a representation in terms of the solution 
of a Painleve V equation |TB|, j3"7f| . 



1 This is a bijection between permutations and pairs (P, Q) of standard Young tableaux with the same 



shape. For expository accounts of the RSK algorithm, see |15, [24| p3| . 

2 Gcsscl [jl6| does not mention random matrices, but in light of well-known formulas in random matrix 
theory relating Toeplitz determinants to expectations over the unitary group, we believe it is fair to say that 
the connection with random matrix theory begins with his discovery. 

3 Ai is the Airy function. For a proof that such a solution exists and is unique, see [|[ [l7| p"3[ . 
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Since the work of Baik, Deift and Johansson, several groups have extended this connection 
between RSK type combinatorics and the distribution functions of random matrix theory. 
The aforementioned result is equivalent to the determination of the limiting distribution of 
the number of boxes in the first row in the RSK correspondence a <-> (P,Q). In the 
same authors show that the limiting distribution of the number of boxes in the second row is 
(when centered and normalized) distributed as the second largest scaled eigenvalue in GUE 
f34jl . They then conjectured that this correspondence extends to all rows. This conjecture 
was recently proved by Okounkov |^8| using topological methods and by Borodin, Okounkov 
and Olshanski @ using analytical/representation theoretic methods. 

Placing restrictions on the permutations a (that they be fixed point free and involu- 
tions), Baik and Rains JO have shown that the limiting laws for the length of the longest 



increasing/decreasing subsequence are now the limiting distributions F\ and F± [36| for the 
scaled largest eigenvalues in the Gaussian Orthogonal Ensemble (GOE) and the Gaussian 
Symplectic Ensemble (GSE). Generalizing to signed permutations and colored permutations 
the present authors and Borodin |J7], showed that the distribution functions of the length 
of the longest increasing subsequence involve the same Fi- Johansson |2l| showed that the 



shape fluctuations of a certain random growth model, again appropriately scaled, converges 
in distribution to F 2 . (This random growth model is intimately related to certain randomly 
growing Young diagrams.) 

Finally, we mention the work of Aldous and Diaconis fTJ where they describe a certain 
one-person card game, called "patience sorting," which is connected to these ideas linking 
Young tableaux with the length of the longest increasing subsequence in either a random 
permutation or a random word; and thence to the limiting distributions of largest eigenvalues. 

At last we come to the subject of the present paper, which is the question of what can be 
said when instead of a random permutation on N letters we have a random word of length 
N from an alphabet of k letters. This may be thought of as a function from {1, 2, • ■ ■ , N} to 
{1, 2, ■ ■ • , k} and it is clear what is meant by a (strictly or weakly) increasing or decreasing 
subsequence. Unlike the case of permutations there is a difference between the two. Given 
such a word w we define ^ N {w) to be the length of the longest weakly increasing subsequence 
in w, and define i^(w) to be the length of the longest strictly decreasing subsequence in 
w.f\ In analogy with permutations we define the distribution functions (giving each word 
probability k~ N ) 

Fj(n; k, N) = Prob fe (4(w) < n) , F D (n; k, N) = Prob fe (e%(w) < n) 
and their generating functions 

oo j.N oo j.N 

G,(n; M) = £ Fi(n; k, N) — G D (n; k, t) = £ F D (n; k, N) — (1.4) 

N=0 iV - N=Q iV - 

Here are our results. We use the standard notation T n (f) for the nxn Toeplitz matrix with 
symbol /. 



4 We could just as well consider weakly increasing subsequences and strictly increasing subsequences since 
reversing the order of the word takes increasing subsequences into decreasing subsequences. In light of the 
RSK algorithm, our choice seems most convenient. 
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Theorem 1. We have 



Gj(n; k, kt) = det T n (fi), G D (n; k, ki) = det T n (f D ), 

where 

fj(z) = /j(z; t) = e*/ 2 (1 + z) fc , / D (z) = f) = e l ' z (1 - 



Theorem 2. G/(n; fc, t) and G?o(n; fe, £) have integral representations in terms of solutions 
of Painleve V equations. 

Theorem 3. Gi(n;k,t) is equal to e fct times the distribution function for the smallest 
eigenvalue in the Laguerre ensemble of k x k matrices associated with the weight function 

Theorem 4. The limiting distribution for the random variable £ T N (w), centered and normal- 
ized, is equal to that for the largest eigenvalue in the Gaussian Unitary Ensemble ensemble 
of k x k hermitian matrices with trace zero.f] More precisely, 

jim Prob, r N{W ) " N/k < s ) = Prob(A 11)nx < ,) . 
1 l2N/k 



The next four sections contain the proofs of these four theorems. 



Theorem 1 is a consequence of Gessel's theorem, just as the permutation analogue (|L1 
is, and the RSK correspondence between words and pairs of Young tableaux. For the con- 
venience of the reader we include a complete proof of Gessel's theorem, containing the main 
ideas of the original but presented somewhat differently. (For related developments, see ||.) 

Theorem 2 is the heart of the paper. The equations are derived very much in the spirit 



of . The logarithmic derivative of the determinant involves a quantity whose derivatives 
in turn involve other quantities. Recursion formulas relating the various quantities allow 
the eventual derivation of a single differential equation which, in the end, turns out to be 
reducible to Painleve V. 

The proof of Theorem 3 consists of showing that the Py function of Theorem 2 for Gi 
is exactly the one which gives the distribution of the smallest eigenvalue in the Laguerre 
ensemble [|35 |. The equation is the same, by inspection, and it is a matter of checking the 



boundary condition at t = 0. 

Given the results of 0] for permutations, it is natural to conjecture that the 

limiting distribution of the number of boxes in the j th row, 2 < j < k, appearing in the 
Young tableaux P in the RSK bijection w <-> (P, Q) is precisely the distribution of the j th 
largest eigenvalue in the finite k x k Hermite ensemble. 

Theorem 4 is proved by an asymptotic evaluation of the multiple integral which gives the 
distribution function for the smallest eigenvalue in the Laguerre ensemble. 

5 The normalization we adopt for the GUE measure is the standard one defined in Mehta p5fl . The 
probability on the right of the displayed formula is the conditional probability given that the matrix from 
GUE has trace zero. 
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After the completion of the original version of this paper there were several relevant 
developments. Johansson 
full conjecture stated above 
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found an independent proof of Theorem 4, in fact of the 
A. Its, using Riemann-Hilbert techniques applied to operator 
equivalents of our Toeplitz matrices, found another proof of Theorem 2. J. L. Snell found 
an error in the original version of Theorem 4. (We thank him for catching the error and 
so saving the authors from further embarrassment.) C. Grinstead found a random walk 
interpretation of the k = 2 problem and used this to determine the limiting distribution in 
this case. And Y. Chen alerted us to the paper [14] of Forrester in which there appeared 
a formula equivalent (given Theorem 1) to the statement of Theorem 3, obtained by using 
identities of Macdonald on hypergeometric functions of several variables. 



II. Gessel's theorem and its specializations 

1. The Cauchy-Binet formula 

We begin by recalling the Cauchy-Binet formula for the determinant of the product of 
two rectangular matrices A and B of sizes mxn and nxm, respectively. We assume n> m. 

Let S mn denote the set of strictly increasing subsequences of length m that can be chosen 
from {1, 2, . . . , n}. For any matrix X of size nxm and any S = {s\, S2, ■ ■ ■ , s m } G S mn , 
denote by X(S\m) the m x m matrix obtained from X by using all m columns of X and 
the m rows numbered by S. Similarly, if X is m x n, denote by X(m\S) the mx m matrix 
obtained from all m rows of X and the columns of X numbered by S. The Cauchy-Binet 
formula is 

det(AB) = det(B(S\m))det(A(m\S)). (2.1) 

ses m n 

Here is a proof. Define 

F(e) =e m det(l + e- l AB) . 
This determinant ismxm and \im. e ^F{e) = det(AB). Then also 

F(e) =e m det (i + e^BA) . 

In the (Fredholm) expansion of this last determinant, the term with coefficient e~ m is 

(BA) ilh (BA) ili2 •■■ (BA) ilim 
(BA) i2h (BA) l2i2 •■■ {BA) l2im 

(BA) imil {BA) imi2 ■■■ {BA)i m i m 

In this sum we can place the restriction that no two indices are equal since when they are the 
determinant is zero. The ml different orderings of . . . , i m } give the same determinant so 
we can drop the ml and sum over all i a with i\ < i 2 < • • • < i m . The terms of order e~i , 
j < m, in the Fredholm expansion do not contribute to the limit as e — > 0. The coefficients 
of the terms e~\ j > m, are zero since the rank of BA is at most m. Finally, each summand 
factors into the product of the two determinants in the Cauchy-Binet formula. 



E 



m\ ■ 

»l,t2,...,tr, 
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The formula remains valid if n = oo if, for example, each row of A and each column of B 
belongs to the sequence space £ 2 . For then AB is well-defined, BA is a finite rank operator 
on £ 2 and the preceding goes through without change. The sum on the right side of ( |2.1| ) 
is then the sum over all increasing subsequences S of length m chosen from the positive 
integers. 



2. Gessel's theorem 



Let V m denote the set of partitions of m, sequences (Ai, A 2 , . . . , A n ) of nonnegative in- 
tegers such that Ai > • • • > A n , Ai + • • • + A n = m, and V = Um>o^m- (Po is the empty 
partition.) We also write A h m when A G V m and denote by £(X) the length of the partition, 
the largest k such that A^ 7^ 0. Let Aq (or A for short) denote the algebra of symmetric func- 
tions over Q. This is a commutative algebra and the vector space direct sum decomposition 
into homogeneous symmetric functions gives Aq the structure of a graded algebra. 

Gessel introduces 

00 

Rn(x,y):= s x{x)s\{y) = £ s *( x ) s ^(y) ( 2 - 2 ) 

A6T> m=0 ASP m 

£(A)<n t(X)<n 

where s\(x) is the Schur function. Gessel's theorem says that R n (x,y) is a Toeplitz deter- 
minant 

R n (x,y) = det(A i . j ) l ^ n (2.3) 

where 

00 

Ai = Ai(x,y) = ^ht +i (x)h e (y) 
1=0 

and h r is the complete symmetric function. (We take h r — for r < 0.) Recall that 

r>0 i>l 



'£hrf = ]l(i-x i ty 



Here is Gessel's proof. Although ( |2.3|) is an identity between two formal power series it 
suffices to prove it when the Xi are real numbers satisfying \xi\ < 1. Let M(x) be the 00 x n 
matrix 

{hi-jix)), (i>l, l<j<n). 

This will be the matrix B of Q2.1|) whereas A will be M t (y). (We interchanged the roles of 
m and n.) Since the entries of the columns of M(x) and rows of M l (y) are exponentially 
small, the formula holds. 

For any increasing subsequence S of positive integers of length n, let Ms(x) be the 
determinant of the n x n minor of M(x) obtained from the rows indexed by the elements of 
S. In the notation of (O) we have M s (x) = det(B(S»), M s (y) = det(A(n\S)). 



Now let A = (Ai, . . . , A n ) be a partition with £(X) < n, and let S = {A„, + i_i + z|l < i < n}. 
Observe that S is an increasing subsequence of postive integers. Then 



M s (x) = det (/iA„ +1 _ 8 - 



-A x , 



l<i,j<n 
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Reversing the order of the rows and columns in this determinant yields 

Ms(x) = det (hxi+j-i(x)) = s\(x), 

where the last equality is the Jacobi-Trudi identity. For such S, summing over all partitions 
A with £(X) < n is the same as summing over all increasing subsequences S of length n. 
Thus 

R n (x,y)=J2M s (x)M s (y), 

s 

and by the Cauchy-Binet formula this is equal to det(M t (y) M(x)). Since 

(M*(y) M(x))ij = E h e-i(v) h t-A x ) = E ht+i-A x ) h e(v) = A i-j, 
i £ 

the theorem follows. 

For the applications which follow it is important to know the the symbol of the Toeplitz 
determinant appearing in (|2.3| ). It is 

oo oo oo oo oo 

<p(z)= E A i (x 1 y)z i = z i J2hi +i (x)h e (y) = Y,he(y) E H^)^ 

i=—oo i=—oo 1=0 £=0 i=—oo 

OO j oo 

= 11 (l-ynz- 1 )' H(l-x n z)- 1 . (2.4) 

n=l n=l 



3. Cauchy's identity from Szego's theorem 

A nice application of Gessel's theorem is a derivation of Cauchy's identity in symmet- 
ric functions^ using the strong Szego's limit theorem for Toeplitz determinants. We have 
(log</?)o = and, for n > 0, 

(logp) n = -E(^) n > 
(io g <fU = -J2(yi) n - 

i>i 

(The subscripts denote Fourier coefficients, as usual.) Applying Szego's theorem (we may 
assume that the X{ and yi are real numbers with absolute value less than 1) then gives 

f 00 1 1 

lim R n (x, y) = exp < ^ - ^ (x { yj) n > = J[ (1 - x^Y > 



n=l i, j>l I i,j 



and hence Cauchy's identity 



3 We freely use various results from symmetric functions which can be found, for example, in J33[, Chp. 7. 
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4. Dual version of Gessel's theorem 



Since ( P73| ) is an identity between two elements of the ring generated by the Xi and yi, any 
endomorphism of this ring yields another identity. Now the complete symmetric functions 
h r are algebraically independent generators of A as are the elementary symmetric functions 
e r . We consider the endomorphism u defined by 

uj(e r ) = h r . 

Then for any partition A = (Ai, A2, . . .) we have 

w(e A ) = h x 

with the usual notation 

e\ = Y[e Xl , h x = '[[hx t . 

i i 

The action on the Schur function is given by oj(s x ) = s x > where A' is the partition conjugate 
to A. 

We define 

Rn{x,y) = u x u y R n (x,y) = Y s \'( x ) s x(y)- 

ASP 
t(\)<n 

Now for any partition A we have £(X) = X[, the length of the first row of the Young diagram 
of shape A', and so 

Rn(x,y)= Y s x'( x ) s \'(y) = Yl s \( x ) s \(y)- 

asp \ep 

(Thus, in this sum we restrict the length of the first row rather than the length of the first 
column.) Applying uj x uj y to (|2.3| ) we obtain the dual version of Gessel's theorem, 

sx(x)s x (y) = det (u x u y A-j(x,y)) i>j=lt _ >n . 

xev 

A j <n 

We record for use below 

00 00 

u x u) y <p(z) = X] e e{y) Yj e i( x ) zt ~ £ 

1=0 i=— 00 

00 00 

= n (1+ynz- 1 ) n(i+x n z). (2.5) 

n=l n=l 

We remark that u y R n (x, y) also equals a Toeplitz determinant. The n — > 00 limit of this 
identity is, by an application of the strong Szego theorem, the dual Cauchy identity 

Sx(x) S y(y) = I]! 1 + XiVj). 
ASP i,j 
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5. Specializations 



If R is a commutative Q-algebra with identity, then a specialization of the ring A is a 
homomorphism ip : A — * R. We always assume that = 1. 

(%) Exponential specialization. If p n = J^a;", the power sum symmetric functions, then the 
exponential specialization is determined by 

ex(p n ) = t5 ln . 

(Recall that the p n 's are algebraically independent generators of A. This homomorphism ex 
is denoted by 9 in Gessel.) 

The fundamental property of this homomorphism is for any symmetric function / 

t n 

ex (f) = ^2[xix 2 - ■ -x n ]f — 

n>0 n - 

where [x±X2 ■ ■ -x n }f denotes the coefficient of X1X2 • ■ -x n in /. Thus if A h N and s\ is the 
Schur function, then 

f N 

«*(»») = /*7n 

where f x is the number of standard Young tableaux of shape A. Hence 



ex 



■ A x2 t 2N 



x ex y R n (x,y) = (f 



N=0 AhJV 
l{X)<n 



(N\) 



2 ' 



By the Robinson-Schensted-Knuth (RSK) bijection 



2 



e(\)<n 



equals the number of permutations <ronJV letters such that £jv(c), the length of the longest 
increasing subsequence in a, is at most n. Hence if each such permutation has probability 
1/iV! we have 

F P (n;N):=Prob(e N (o)<n) = ±- ] £ (f x f . 



AhJV 
£{\)<n 



Thus we know that its generating function is given by 

oo ^2N 

F P (n;N) — = ex x ex y R n (x,y). 
n=o iV ■ 

Gessel's theorem tells us that R n (x,y) is the n x n Toeplitz determinant with symbol 
(p(z) given by (|2.4j ). This may be written 

<P( Z ) = J2 h r(y)z~ r J2 h s(x)z s . 

r>0 s>0 
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The important observation is that since ex is a homomorphism, ex x ex y R n (x, y) is the 
Toeplitz determinant with symbol 

f P (z) := ex x ex y { V {z)) = £ £ = e'^l 

This is precisely ( |1 . 1| ) after changing t to \fi. 

(ii) Principal specializations. The principal specialization of order n of g is defined by 

PSn(g) = g(l,q,q 2 ,...,q n -\0,0,...). 

(Thus we replace Xi by q l if i < n and by itherwise. If we let n —>■ oo we obtain the stable 
principal specialization of /.) Setting q = 1 in ps n gives 

ps 1 n ( J g)= g(l,l,..., 1,0,0,...) 

where 1 appears n times. 
Observe that 

p4( s a) = 4(*0 

where d\(k) is the number of semistandard Young tableaux of shape A that can be formed 
from an alphabet of k letters. (Recall that a semistandard tableau is weakly increasing across 
rows and strictly increasing down columns; a standard tableau is strictly increasing across 
rows.) This is most easily seen from the combinatorial definition of the Schur functions. 
Applying the homomorphisms ps\ and ex to R n gives 

oo / 

(psl) x ex y R n (x, y) = £ ]T d x (k) f x I ' (2.6) 



N=0 \ AhJV 
\e(\)<n 




The RSK correspondence associates to each word w of length iV formed from an alphabet of 
k letters a pair of tableaux, (P,Q). Here the P are semistandard Young tableaux of shape 
X\- N made from the alphabet {1,2,..., A:} and the Q are standard Young tableaux of shape 
A h iV on the numbers {1,2,..., iV}. Thus 

XhN 
e(\)<n 

counts the number of words w of length iV with strictly decreasing subsequences all of length 
less than or equal to n. Obviously, 

E d ^ k ) f x = kN - 

XhN 

The symbol of the Toeplitz determinant that equals the generating function (|2.6| ) is 

k £r e t/z 

Sd{z) = Id{z; t) := (ps{) x ex y ((p(z)) = JJ(1 ~ z )~ l E ~\ z ~ r = n _ 

i=l r>0 r - V 1 Z > 
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Hence we have shown that if £n( w ) denotes the length of the longest strictly decreasing 
subsequence in word w, and if each word of length N is assigned probability 1/k , then the 
generating function of the distribution function 

F D {n; k, N) := Prob fe (f%(w) < n) 

is given by the Toeplitz determinant with symbol fo'- 

Y: F D (n; k, N) = det (T n (f D (z; t)) . 

N=0 iV • 



Recalling ( |1.4| ) we see that this is 

G D (n;k,kt) = det (T n (f D (z;t)). (2.7) 

Since, under general conditions, changing the symbol of a Toeplitz matrix from f(z) to 
f(az) is a similarity transformation, the associated Toeplitz determinant does not change. 
Therefore the symbol of the Toeplitz determinant in ( |2.7p may be replaced by 



pVt/z 

f(Viz/k; t/k) 



;i - Vtz/k) k 

whose k — > oo limit is e^^ z+z \ This shows that for fixed N, 

lim F D (n; k, N) = F P (n; N). 



k^oo 



Again, this is intuitively clear since as the size of the alphabet approaches infinity, any 
random word of length iV is very likely a permutation. (This also uses the fact that the 
distribution of the length of the longest decreasing subsequence of a random permutation is 
the same as the distribution of the length of the longest increasing subsequence.) 

Finally we apply the same specialization {ps\) x ex y to the dual version of Gessel's theo- 
rem. We see that ( |2.6| ) is replaced by 

00 I \ t N 

(psl) x ex y R n (x,y) = £ I E d x (k)f x \ — . (2.8) 



N=0 \ Aj<n 

xev 



(We used here the fact that f x ' = f x , which follows immediately from the hook length 
formula for / A .) Thus we obtain the generating function for the distribution of the length 
of the longest weakly increasing subsequence. Using ( |2.5|) we find that the symbol of the 
Toeplitz determinant that gives ( |2.8|) is 



fi(z) = fi(z; t) := (psl) x ex y J2 e r(x)z r £e s (y) 



r/ z 

r>0 s>0 



(psl) x ex x J] (1 + xjz) £ e s {y)z s = (1 + zf e t/z . 

j s>0 
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To summarize, we have shown that if £^{vj) denotes the length of the longest weakly 
increasing subsequence in word w of length N, and if each such word has probability l/k N , 
then the generating function of the distribution function 

F/(n; k,N) := Prob fe (l»<n) 

is given by the Toeplitz determinant with symbol //. Precisely, 



Gi{n;k,kt) = Fj(n;k,N) 



(kt) 



N 



N=0 



det (T n (fi(z; t)) . 



The same k — > oo remarks hold here as in the strictly decreasing case. 
Relations ( |2.7[ ) and ( |2.9| ) are the assertions of Theorem 1. 



(2.9) 



III. Recursion and differentiation formulas 

1. Universal recursion relations 

In this section / will be an arbitrary function, with Fourier coefficients f\ and associated 
n x n Toeplitz matrix 

T n (f) = (fi-j), (i,j = 0,---,n-l). 

We assume T n (f) is invertible and obtain several relations connecting various inner products 
involving T n (f)~ l . Most of these relations actually appeared in P7|. There our T n (f) was 
symmetric and unfortunately some of the relations derived in [|37| used this fact. Since this 
does not happen here, everything has to be modified for the more general case. A reader of 
the earlier article will find familiar much of what we now do. 

We introduce the n-vectors 



/ 1 \ 






/ \ 




v 1 ; 



r 



t h \ 

h 

fn-l 
V fn J 



f- 



( fn \ 

fn-l 

h 

V h J 



and define / by f(z) — f(z 1 ), so that T n (f) is the transpose of T n (f). We write 

A = T n (z- 1 ), A' = T n (z). 

Thus A is the backward shift and A' is the forward shift. It is easy to see that 

T n {z- X f) = T n (f) A + /+ ® 5 + = AT n (f) + r ® r , 

T n (z f) = T n {f) A' + f-®6- = A' T n (f) + 6 + ® /+. 
These identities explain why the vectors / and / arise. 



(3.1) 
(3.2) 
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The inner products involving T n (f) 1 are 

If one of these quantities defined in terms of / is given a tilde, then / is to be replaced by / 
everywhere in its definition. Thus, for example, 

Note that V+ = V+ = D n /D n+1 , where 

D n = D n (f)=detT n (f)- 

Some other inner products may be expressed in terms of these using the isometry that 
reverses the order of the components of a vector (and replaces a Toeplitz matrix by its 
transpose). Thus, for example, 



(TnC/)" 1 /-, 5 + ) = (T n (f)- 1 f+, 5-) = U- 



(3.3) 



We shall use this isometry from time to time below without comment. 

The basis for all the universal relations we shall obtain is the following formula for the 
inverse of a 2 x 2 block matrix: 



/ A B 



\C D 



(A — BD~ 1 C)~ 1 



(3.4) 



Here we assume A and D are square and the various inverses exist. Only one block of the 
inverse is displayed and the formula shows that A — BD~ 1 C equals the inverse of this block 
of the inverse matrix. 

We apply fl3.4|) first to the (n + 1) x (n + 1) matrix 



/ 

fi fo 

V fn fn—1 



1 \ 

/-n+1 



with A = (0), D = T n (f), B = (0 ••• 1), C = f+. In this case A - BD~ l C = 
~{T n {f)~ 1 f + , 8~) = —U~. This equals the reciprocal of the upper-left entry of the inverse 
matrix, which in turn equals (— l) n times the lower-left n x n subdeterminant divided by 



D n . Replacing the first row by (fo /_] 



/ fo f-i 
fi fo 



V fn fn—1 



f- n ) gives the matrix 



f-n \ 

/-n+1 



T n+ i(f). 



(3.5) 
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The lower-left entry of its inverse equals on the one hand (T n (f)~ 1 5 + , 8~) = V~ +1 and on 
the other hand (— l) n times the same sub determinant as arose above divided by U~ +l . This 
gives the identity 

D n +i V n +i 



u- = K+i 



D r , 



(3.6) 



If we take A to be the upper- left corner of ( |3.5|) and D the complementary T n (f) then 
C = f + and B = ■ ■ ■ /_ n ), and we deduce that 



fo - (T n (/)-V + , /+) 



(3.7) 



n+l 



Next we consider 



/ /q /-i 

/i /o 

fn fn—l 
\ fn+1 fn 



f—n f—n—1 ^ 

/—n+l jf— n 

fo f-1 

fl fo J 



T n+ 2{f). 



We apply to this an obvious modification of ( |3.4j ), where A is the 2x2 matrix consisting of 
the four corners of the large matrix, D is the central T n (f), C consists of the two columns 
/ + and /~ and B consists of the rows which are the transposes of / + and Then we find 



A-BD^C 



( fo - (Tn(f)- l f + , f + ) f-n-l ~ (Tnif)- 1 ?-, /+) 



V - (r n (/)-v + , n fo - (T n (fr i r, /-) 

and our formula tells us that this is the inverse of 



V KT+2 v T 



n+2 



This gives the two formulas 



fo-(T n (f)- l f + ,f 



V, 



n+2 



Ki+2 Ki+2^n+2 



fn+1 — (T n (f) X f + , f 



n+2 



Vn+2 Vn+2^n+2 



Comparing the first with (|3.7|) we see that 



V + —V~V~ — v + v + 

v n+2 v n+2 v n+2 ~ v n+l v n+2i 

and therefore that the preceding relations can be written 

1 



fo ~ (T n (/)-V + , / + ) 



V+ ' 

v n+l 



fn+i - (T n (/)- i / + , n 



1 v. 



n+2 



Vn+l Vn+2 



(3.8) 



(3.9) 
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Notice that Q) and give 



(3.10) 



n+l 



Next, we apply (|3.4j ) to the matrix 



/ 1 

fi fo 

V fn fn 



o \ 

7-71+1 

/o J 



Now we have A = (0), D = T n (f), B = (1 • • • 0), C = f + and so A - RD^C = 
— (T n (/) _1 / + , 5 + ) = — Therefore this equals the inverse of the 0, entry of the inverse 
of the matrix, which in turn equals its determinant divided by D n . But its determinant 
equals —K, where K is the 0, 1 cofactor. Thus 



^n 



Now look at the matrix ( |3.5| ) and consider the 1, entry of its inverse. It equals on the one 
hand (T n+ i(/)) _1 5 + , A'5 + ) and on the other hand —K/D n+ i. This gives the identity 



u: 



D 



n+l 



(T n+1 (f)- 1 S + ,A'S + ) 



Dn 



V, 



(T n+1 (f)- 1 6 + ,A'5 + ). 



(3.11) 



n+l 



To evaluate the inner product on the right side we for the moment replace n + 1 by n so 
that we can apply earlier formulas. The inner product becomes 

(T n (f)~ 1 5 + , A'5 + ) = (AT n (f)~ 1 5 + , 5 + ). 

Multiplying the second identity of ( p.l|) left and right by T^/)" 1 gives 

A TrXfY 1 + T n {fY l f + ® T n {fy l 5 + = TM- 1 A + T n (f)~ 1 5~ ® T^f)- 1 /-. (3.12) 

Applying this to 5 + (observing that A5 + = 0) and taking the inner product with 5 + we find 
that 

(A T n (fy 1 5 + , 5 + ) = V- U- - U+ V+. 
Here we used (|3~3|). Replacing n by n + 1, we see that ( |3.11| ) becomes 



V. 



n+l 



u: 



y/+ w n+l 
^n+l 



u, 



n+l- 



This gives 



by (ED- 



U n U n+1 — U n U n+1 



(3.13) 
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Next we apply ( |3.4| ) to the n x n matrix 

(0 f-i 
f Q 



f-n \ 
f-n+1 



f-1 

fo J 



1 fn-2 

\0 / n _i 

The formula says that the inverse of the upper-left entry of its inverse equals 

-(T n (f)- 1 A8-J + ). 



This inverse also equals (— l) n 1 D n /K where K is the n — 1,0 cof actor of the matrix ( |3.5| ). 
But (T n+1 (f)~ 1 A5~, S + ), the 0, n - 1 entry of the inverse of (5), equals (-l) n ~ l K / D n+l . 
We have shown 

~ (T n+1 (f)- l A5-, 5 + ) = -(T n (f)- l A5-, /+). 



"Or 



This may be written 

(5 + , AT n+1 (/) _1 5~) = -V n+1 (5 + , AT^f)- 1 ]-). 
To compute the left side we use the fact that 

/ /-n+l \ 

Af = 



f-i 

V o j 



T n (f) S~-fo6- 



Hence 



TM^Af- = S--f T n (f)- 1 6-. 
Therefore applying ( |3.12| ) to /~ gives 

AT n {f)- 1 j- = 5--hT n {f)-H- 

+(T n (/)- 1 r, I') T n tf)-H- - (T n (f)- l 5 + , /-) T„(/)-V + . 
Taking inner products with 5 + we find 

(6 + , atm)- 1 n = - (/o - (T n (f)- i r, r) v- - u- u+. 

But Q3.9p gives 

1 



V. 



n+l 



Hence 



AT n (/) -1 /" 



w n w n ' 



n+l 



(S+, AT n+1 (f)- 1 5-) = V- + U~ U+ V+ +1 . 



(3.14) 
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2. Nonuniversal recursion relations 



Here we restrict to our symbol 

friz) = e*/* (1 + z) k 

but write / instead of fi for notational convenience. We shall obtain relations which follow 
from the representation 

2m J 

upon integrating by parts. The fact 

— [ ^-\e^(l + z) k+1 z^} dz = 
2mJdz^ v ; J 

gives 

/ (f j+ i + /;) + (k + l) ./• . - j (./• + ./• ,) = 0, 
U + t)fj + U-k-l)fj-i+tf j+ i = 0. 
Replacing j by i — j we obtain the i,j entry of a matrix identity: 

(M + t) T n (f) - T n (f) M + (M - k - 1) T n (*/) - T n (z/) M + tT(z~ 1 f) = 0. 

Here M denotes the diagonal matrix with diagonal entries 1, 2, • • • , n. We use the identities 



( |3.1|) and (|3.2| ) to write this as 

(M + t) T n (f) - T n (f) M + (M - k - 1) (A' r n (/) + 5+ <g> /+) 

-(T n (/) A' + r ® r)M + t (T n (f) A + /+ ® 5+) = 0, 
(M + t) r n (/) - r n (/) M + (1 - fc) A' T n (/) - T n (/) A' + tT n (f) A 
-A; 5+ ® f + - n f~ ® <T + t /+ ® 5+ = 0. 
We multiply this left and right by T n (f) , obtaining 

TM)- 1 (M + t) - MT n (f)~ l + TMr 1 {M — k — l)M — A'MTnif)- 1 +tAT n (f)- 1 

-A;T„(/)- 1 ( 5 + ®T„(/l-7 + -nT„(/)-7-®T„(/)- 1 r+tT„(/)- 1 / + ®T„,(/)- 1 5 + = 0. (3.15) 

This is the basic matrix identity. Applying this matrix identity to 8 and taking inner 
products with 5 ± gives identities for scalar quantities. We shall need three of the four. 

First we apply ( |3.15|) to 5 + and take the inner product with 5 + . If we recall our definitions 
and the fact A5 + = we obtain 

t V n + + (1 - k) (Uf)- 1 A' 5 + , 5 + ) + t (A T n (f)~ 1 8 + , 5+) 

- k V+ U+ - n (T^/)" 1 /- , 5 + ) V- + t U+V+ = 0. (3.16) 

But 

[T n {f)- i r, 5 + ) = (T n (/)-v + , n = u-, 
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so ( |3.6| ) gives 
Thus 



v: 



-v n + u-_ v 



(TMT'r, <*+) v~ = -v n + u- u-_ v 

Next, by ( |3.11| ) we have 

(T n (f)- 1 A'5 + , 5 + ) = -V+U^, (AT n (fy 1 5 + , 5+) = -V+U+.V 

Substituting these relations into ( |3.16|) and dividing by gives 

t + (k - 1) U+_ x - t U+_, -kU++n U- U~_ x +tU+ = 0. 

But ( |3.13| ) gives U~ U~_ l = U£_x — U£, from which we see that the above becomes 

t + (k + n-l)i^ 1 -(k + n )i^ + t(U+-U^ 1 ) = 0. 

The derivation holds also for n = 1 when one defines Uq = Uq = 0, as one can easily check. 
Therefore summing over n gives 

nt- (k + n)U+ + tU+ = 0. (3.17) 

For the next relation we apply (|3.15 ) to 5~ and take the inner product with 5 + to obtain 

(n + 1 - 1) V- + 1 (A Uf)- 1 5-, 5+) - k V+ U- - n U~ V+ + t U+ V~ = 0. 

If we applying ( |3.14| ), divide by and use ( |3.6| ) we obtain 



-{n + t-l)U-_ x + t 



v- 



+ UU U-_, \-(k + n)U--t U+ U~_, = 0. 



But 



n-l 



~U n -2 (1 U n _ l U n _ l/ 



V + V + i v + 

n r n—1 r n 

by ( p.6| ) and ( |3.1(J| ). Replacing n by n + 1 and introducing the function 



^ n = i-u~u- 



we obtain 



(t + n) U~ + t ([/+ U~ - U-_, $ n ) - (k + n + 1) U~ +1 - t U n 



0. 



By fl3.13p this may be written 

- (t + n) U~ + t {U-f U- +l - t U~_ x <S> n -(k + n+ 1) U' +1 = 0. 

Finally we we apply fl3.15| ) to 5~ and take the inner product with 5~ to obtain 
n V n + -(n-l) (T n (f)- l 6-, A5~) — k V~ U~ - n U+ V+ + 1 U~ V~ = 0. 



(3.18) 



From ( |3.11| ) we see that (T n (f) 1 5 , A 5 ) = -f/^„i V n + . Substituting this into the above 
identity, dividing by V n and using (|3.6| ) gives 

t + (n - 1) U+_ x + fc Ed E/- -nU+-t U~ fj-_ t = 0. (3.19) 

Using (|3.13| ) we can write 

k U-_, U- = n U-_ x U- + (k-n) (U^-i ~ 

Substituting this into the preceding and replacing n by n + 1 give 

t + (k — l) U+ + (n + 1) C/- - k U+ +1 - t U~ +1 U~ = 0. (3.20) 



Another identity can be gotten by using (|3.13|) to replace n {U n _i — U n ), which we see 
in ( |3.19|) , by nU^-i U~ . This gives 



* " Uti + (* + n) U~_ x U- - t U- U-_, = 0. 
Replacing n by n + 1, multiplying by C/~ and adding to ( |3.1S| ) gives 

nU- + U-U+ + $ n ((k + n + I) U~ +1 +t [/-_,) = 0. 
3. Differentiation formulas 



(3.21) 



We continue to take f(z) = e^ z (1 + z) k and write D n (t) for D n (f) Since d/*/c?t = z 1 f 
we have 

^ log D n (*) = tr Tn(f)- 1 T^z- 1 /) = tr (A + T^f)' 1 f + ® 5 + " 
by and so 



d 
It 



log D n (t) = W n 



(3.22) 



which is why this quantity arises. Others will arise from further differentiation. 
We use 



d_ 

dt 



Tn(f)- 1 = -T n (/)- 1 r n (z- 1 /)T n (/)- 1 = -AT„(/) _1 — T n (f)~ 1 f + <g> T ri (f)~ 1 5 + 



df+ 
dt 



[z- 1 f)+ = Af+ + f n+1 5-. 



Hence 



dU+ 



= -(AT n (f)- l f+, 6 + ) - (T n (f)-'f + , 5+) (T n (f)-'S + , /+) 



-i ^+ 



dt 



+(T„(/)- 1 A /+, 5+) + / n+1 (TM^S-, 5 + ). 
By ( |3.12|) the two terms involving A combine to give 

(T n (/)-V + , 5+) (T n (/)-M + , /+) - (T n (f)~~ 1 8~ , 5+) (T n (f)-if-, f 
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Using this we find that the preceding simplifies to 



-\C(T n (/)- 1 r,/ + ) + /n+iV; 



-v- J- Y&i 

n V+ V + ' 

v n+l v n+2 



by the second part of ( |3.9| ). This equals 



r n r n v n+2 



V+ V"tLi U H 



+2 



by ( p.6|) and ( [3.10] ). We have shown 

dt 



(3.23) 



In completely analogous fashion (we spare the reader the details) we compute 

<Hr -v n + (T n (f)- i r,f + )+f n+1 v+, 



dt 



and using again the second part of (|3.9p, ( |3.6| ) and ( [3.10[ ) we find that 

dU: 



dt 



(3.24) 



To find formulas for the derivatives of we use 



dt 



Tnifr 1 = -T n (/) -1 T n (zf) T„(/) -1 = — A'T n (/) _1 + T n (f)~ l f~ ® T n (fy l 5~ 



df+ 
dt 



(zf) + = A'f + + f S + . 



At the appropriate points in the computations we use the analogue of (|3.12| ) with A' instead 
of A, and the first part of (|3.9| ) rather than the second. Again we spare the reader the details. 
The results are somewhat simpler: 

dU+ 



- $ 

dt n ' 



dt 



-& n U n _ v 



(3.25) 
(3.26) 



Observe that from ( |3~22D , Q3T7|) , and Q3T23|) we have 

j^t j log D n (t) = (k + n) $ n - n. 
This gives the representation 

log D n (t) = (k + n) [ log(t/t') $„(t') eft' - n*. 
jo 
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IV. Painleve V and the Laguerre ensemble 



1. Derivation of the differential equation 



We begin by showing how differentiation formulas ( |3.24|) and ( ]3.26|) have analogues in 
which only indices n and n — 1 appear on the right side of ( |3.24j ) and only indices n and 
n+1 appear on the right side of (|3.26|) . 

Solve (|3.17| ) for U~ +l . The solution involves U^+i, which we can solve for in ( |3.20| ). Thus 
U~ +l , and so dU~ /dt, can be expressed in terms of quantities with indices n or n — 1. To 
obtain a differentiation formula for U~ that involves only n and n + 1 simply solve ( |3.20|) 
for Un-i- The results of this are 



dt 



-u- 

t n 



t 



: I 12— TJ- TJ~ 2 



dt 



n ~ 1 ~ ~ 1 

-U n +-U+U- + ~(k+l + n)$ n U- +1 . 



Now compute $'„ using (|3.24 ) and (|4.2j ). The result can be put in the form 
K n* n -l l U+{* n -l) _ u-- l -(n + k + l)U- U~ 



(4.1) 
(4.2) 

(4.3) 



Now solve ( |3.18| ) for U n+1 U n and insert the result in the right hand side of ( |4.3| ). The result 
can be written as 



kU 1 t +1 -(2n + 2 + k)U;U- +1 =t 



n 



n 



$ n -l 17+(S„-1) 



t-u: = s, 



say. Noting that 
we have 



u n u n+l ~ u n u n+1 



u- u- +1/ 

by ( p.25|) . Therefore differentiating the left side of (|4.4|) gives 

£' = -(2n + 2 + £;)$„, + 2(& + n + 1) $ n+ i, 

2(fc + n + 1) $ n+ i = £ ' + (2n + 2 + k) $ n . 
Computing £' using the right side of Q4.4J) we obtain the representation 



n j. V n 



/ 



2(n + k + 1) = 2 + 2(n + fc) - (n + ttf - $ n )-f - t 



$2 



Simply integrating the preceding equation using ( |3.25|) gives 



2(n + fc + i)^ 1 = t|^-n%^ + 2(fc + n)^ + ^ + t. 



3: r. 



(4.4) 
(4.5) 



(4.6) 



(4.7) 
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We now write out the last relation from which the differential equation will follow. 

(u~ +1 u-) (u-u- +1 ) _ 1 (u~ +1 u-) u-u- +1 ) 



$ n+1 = 1 - u n+1 u n+1 = 1 



u-u- 



l-$ r 



Use ( |4.5| ) and (|4.7|) to express U~ U^+i * n terms of $ n and $' n . Similarly use ( |4.3[ ) 
and (|4.7| ) to express C/~ in terms of the same quantities. Finally, on the left hand side 
of ( |4.8| ) use (f4.6|). The result is a third-order differential equation for [/"+. (Recall ( |3.25| ).) 
This third-order equation is 



1 



w 



2 \w' w' - 1 
t + n 



w 



n\2 



W 



n-t + 2w) - 



t t 
(n + w) 2 



2(k + n) , 
■ w 



t ' ' 2t 2 v ' 2t 2 w> 

Cosgrove tells usQ that the equation integrates to 



2t 2 (w' - 



(4.9) 



t 2 ( w "f = -A(k + n) t (w'f + [A(k + n)w + t 2 + 2(2k + 3n) t + n 2 } (w 1 ) 
- [2(t + 2k + 3n) w + 2nt + 2n 2 } w' + (w + n) 2 . 



(4.10) 



The a = a(t) form of Painleve V as given by equation (C.45) in Jimbo-Miwa [|19| (see 
also p7| ) is, after changing a to —a and taking the special parameter values v$ = v x = 0, 
^2 = k, u 3 = k + n, 



(t a"f = [a - t a' - 2 (a') 2 + {2k + n) a'} 2 - 4 (a') 2 {a' -k)(a'-k-n). 



(4.11) 



If 



w = t 



a 



(k + n) 



then ( |4.10| ) and (|4.11| ) are equivalent. Notice that since w 
kt — tU+ and therefore by QX22; ) 

d 



and therefore 



(|3.17|) says that a 



a = -t-\og (e~ kt D ri 
e~ kt D n (t) = exp (-jl^T-dA 



(4.12) 



This, with Theorem 1, gives Theorem 2 for Gj. For Go it is simply a matter of changing 
k to —k and t to — t. 



7 Cosgrove, in his analysis of certain third-order differential equations, has shown that the third order 
differential equation of Chazy Class I (see (A. 3) in |1(|) can be integrated to a second order and second 
degree "master Painleve equation" (see (A. 21) in [jl0|). This master Painleve equation, called SD-1 in ||], 
contains all the Painleve equations I- VI. O ur (|4.S| ) is a special case of Cosgrove's (A. 3). Carrying out this 
reduction [[H] in this special case results in (4.10). Cosgrove's integration constant equals — n 2 /4 in our case. 
This follows from the boundary conditions derived below in (4.14). 
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2. Laguerre ensemble interpretation of D n (t) 



In order to specify which solution of (|4.11 ) our a is, we must determine the boundary 



condition a satisfies at t — 0. We have by (|3.25|) and ( p^6| 



Now is the upper-left entry of T„ +1 (/) _1 where, recall, f(z) = e t//z (1 + z) k . As t — > 
this approaches the upper-left entry of (/ + A')~ fc , which is clearly equal to 1. Equally clearly, 
the lower-left entry of the inverse has limit (~~^) , so that 

To determine the behavior of V~ +1 , the the upper-right entry of T n+ i(/) _1 , we write 

T n+1 (f) = (I + A') k + £ ^(*W«>, (4.13) 

p>0,q>0 P- \1J 

where A^ denotes A J if j < and (A') - - 7 if j < 0. Factoring out (/ + A') k and taking the 
inverse gives 

-l 

-1 _ f r i ( t i a '\— Jfe ^ / ^ \ a (p-o) I / r i a 



T^if)- 1 = \ I + (I + AT k E Zt(J A(P ^I ( J + A 

p>o 

If we expand out the inverses we get a sum of products. Each product has factors of the 
form t Pi A( pi ~ qi ^ and other factors which are nonnegative powers of A'. Such a product will 
have a nonzero upper-right entry only if X)(Pi + Qi) > n- Therefore, since each pi > 1, the 
lowest power of t which can occur is n. Moreover this power occurs only when all qi = 
and all the nonegative powers of A' which occur in the product are 0. This means that we 
get the same lowest power of t in the upper-right entry of the inverse if in (}4.13|) we replace 



(/ + A') fc by / and in the sum we only take the terms with q = 0. This amounts to replacing 

W/)by 



Eh* 



e 



t A 



p>0 



pi 



The inverse of this operator is e tA and the upper- right corner of this matrix is exactly 
(-l) n t n /n\. Thus 



(—t) n 

v- +1 = [ —r + o(t n+l : 



as t — > 0, and so 



% = (k + n) (~ k ) + (t^) = ^( n + h )t n + O(f^) 

at \ n n\ n! V n 
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Since <r(0) = 0, 

Here is the remarkable fact: the same function a which satisfies the equation (|4.11 ) to- 
gether with the boundary condition ( [4. 14] ) gives a representation for the Fredholm determinat 
which equals the distribution function for the smallest eigenvalue in the Laguerre ensemble 
of k x k matrices associated with the weight function x n e~ x . Precisely, we have 

Prob (A min >t) = det (J - K L ) , 

where is the integral operator on (0, t) with kernel 

K L (x, y) = [k(k + n)] 1 ' 2 Vr>A*)VL,k-i(v)-VL,k-i(x)<p L , k (y)^ 

x-y 

Here 



Moreover 



det (/ - K L ) = exp (- jf* ^p- dt'j . (4.15) 

(See [[55[], Section VB.) It follows from this and (|4.12| ) that 

e- kt D n (t) = det {I-K L ), 
which, with Theorem 1, is the assertion of Theorem 3. 

3. Limiting distribution as iV — > oo 

The foregoing can be restated in more concrete terms as 

Fi{n\ k, N) = e kt c Kn [°° ■ ■ ■ [°° U x] e~2> A(x) 2 dx x --- dx k , (4.16) 

where A(x) = Yii<j{ x i ~ x j) an d c k<n is the normalization constant defined by 

/■CO roo ^ 

( c *,0 =/ ■'■J l[x]e-^^A(x) 2 dx 1 ---dx k . 
In fact (g§, formula (17.6.5)) 

(c fc , n )- 1 = l!2!---fc! Y[(n + j)\. (4.17) 

3=0 

If we make the variable changes Xj — > Xj + t in the integral, the right side of fl4.16|) 

beCOmeS ,00 roo 

c k ,n ■•■/ Y[(xj + t) n e A(xfdx l ---dx k . 
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Therefore 

F I (n;k,N) = - j ^c k , n J q e-^Aixydx—jt-^lKxj + t^dt, 

where the inner integral is taken over a contour surrounding t — and we write <ix for 

C&Ti • • • C&Tfc. 

Set 

N = kn-r, r = [skV2n]. (4.18) 

Then 

± J H( Xj +tTdt=±J f exp {n £ log (1 + r 1 x 3 )} * 

= i / exp { r logt + n ( r 1 S *' " 2^ 2 E *? + •••)} * 

as long as the integration is over a contour where \t\ > J2 x j- 111 fact we integrate over the 
circle \t\ = n J2 x j/ r - The function rlogt + nt~ x J2 x j has a critical point at t — n J2 x j/ r 
and its real part restricted to the circle has an absolute maximum there. The rest of the 
exponent, n y—\t~ 2 Y, x ^ + • • •), is uniformly bounded on the circle and equals 

at the critical point. It follows that as n — > oo we have, uniformly for all Xj, 



r 2 E^- 1 (nE4 



eXP \ 2n (E^) 2 J r(r + l) • 
Thus 

*. *0 ~ ^ / • • •/«" {-s J^} E*i)'«- s,J A <*> 2 

(4.19) 

Define 

Z:={(x,)ER k : £x,=0} 
and for general (xj) G R fc write 

2/ = 5Z^i> xj = -x'j+y/k, 

so that (x^) G We integrate over Z with Lebesgue measure and over y G R. Since each 
> the y integration is restricted to 

y > k maxx'. 
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We find that the double integral in ( 4.19[) equals (after changing back from x' to x) 



■Jk. ' Z Jk maxrr, I 2n y 

r 2 /2kn . ( f 2 J- X 2 ) 

A(x) 2 dx I exp < — ^- > y r e~ y dy 



y/k J Z Jk max Xj 1 2?2 y^ , 

-r 2 /2kn r roo ( 2 y T 2 

= ^_(2n)(^)/ 2 / A(*) 2 d* / exp -f , ^ 

V/C ^ « maxij [ An y A J 

where to obtain the last we made the substitutions Xj — > y/2nxj, y — > V2ny. 

The factor y r e _y in the inner integral achieves its maximum on R + when y = 
r/\^2n = sk + o(l), at which point the other factor in the integral equals e - ^^. Hence if 
maxxj < s (so that sk is interior to the range of the y integration) the inner intregral is 
asymptotically equal to 

e -£*| I™ tfe-v^dy = e"M (2n)- (r+1)/2 r(r + 1), 
Jo 

while if maxxj > s the inner integral is o of this. Moreover the inner integral is at most 
/ y r e- yV ^dy < e " fcmax ^ / y r e ^^-i) dy < c e -k max., ^n)- {r+1)/2 T(r + 1) 

J fc max a; j J 

for a constant C independent of the Xj and n. Hence application of the dominated conver- 
gence theorem shows that if we define 

Z s = {x G Z : maxxj < s} 



then the double integral in (|4.19|) is asymptotically 

-r 2 /2kn 



( 2n )( fc2 - 1 )/ 2 r(r + l) / A(x) 2 dx. 

JZs 



If we recall the definition ( |4. 18| ) of r and the value of Cfc n given by ( |4 . 1 7] ) and apply Stirling's 
theorem we obtain 

FAn; k, N) ~ 7fc / e~5>; A(x) 2 dx, 
Jz s 

where 

% 1 = l\2\---k\ (2nY k - 1 ^ 2 2-( fc2 ~ 1 )/ 2 . (4.20) 

Equivalently, 



ri , ,£ I N (w)-N/k 
lim Prob fe I ^ v / '— < s 



2N/k 



j = 7fc jf A(x) 2 cfe. (4.21) 



The right side is the conditional probability that the largest eigenvalue of a matrix from 
GUE is at most s, given that the matrix has trace zero, and so Theorem 4 is established. 
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4. Large k asymptotics 



In this section we denote by F(s, k) the probability in k x k GUE that A max < s and by 
F°{s,k) the same probability in k x k traceless GUE. Thus F°(s, k) is given by the right 
side of Q4.2 H) while 

F(s,k)=c k [ ■■•( e~S^ 2 A(x) 2 cfei ■ ■ -dx k , 

J—oo J—oo 

where (see ( 14.20 ) and formula (17.6.7) of ||25|| ) c k = 'Jk/V^- Using now the variable change 

2/ = XX' x j =x' j + y/k 
and then changing x' back to x as before we find that 



F(s, k) 



Ik 
wnk 



dy 



-S-f+V) A(x) 2 dx 



-y/k 



r e'^dy [ e-^ x >A(x) 2 dx=^ f°° e"^ 2 F° 

IX J-oo JZ s -y V 7f J-oo 



(s - y, k) dy. 

Replacing s by \/2k + s/\/2k 1 ^ and making the variable change y — > yj\p2k v ^ we obtain 

fc l/3 



F^ + s/V^A; 1 / 6 , fc) = — p= /°° e-^Pl^iVU+is-y)/^' 6 , k)dy. 

\/27T J-oo 



The factor F°(\/2k + (s — y) / \/2k 1 ^ , k) in the integrand is a bounded and nonincreasing 
function of y whereas the integral over \y\ > 5 is o(l) for any 5 and the right side without 
this factor equals one. Hence for any 5 > 

F°{V2k + (s- S)/V2k 1/6 , k) + oil) < F(V2k + s/v^k 1 ^, k) 

< F\V2k + {s + 5)/V2k 1/6 , k) + oil). 
It follows from this, (FT73) and the fact that F 2 (s) is continuous that 

lim F°{V2k + s/V2k 1/6 , k) = F 2 {s). (4.22) 

fc^oo 



(A result which includes this, where N — ► oo and k — > oo with N ^> k, can be found in |22 



Another way of expressing ( 4.22|) is as follows. Let i k equal the weak limit of 
(£jy — N/k)/y/N as N — > oo. The distribution function of is 

lim Prob fc JVV 7 < s 



F"(sJk/2, k). 



Then fl4.22|) is equivalent to the statement 

lim Prob ((4 - 2)0 3 < s) = F 2 (s). 
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